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At non-equilibrium phase transitions into absorbing (trapped) states, it is well known that the 
directed percolation (DP) critical scaling is shared by two classes of models with a single (S) absorb- 
ing state and with infinitely many (IM) absorbing states. We study the crossover behavior in one 
dimension, arising from a considerable reduction of the number of absorbing states (typically from 
the IM-type to the S-type DP models), by following two different (excitatory or inhibitory) routes 
which make the auxiliary field density abruptly jump at the crossover. Along the excitatory route, 
the system becomes overly activated even for an infinitesimal perturbation and its crossover becomes 
discontinuous. Along the inhibitory route, we find continuous crossover with the universal crossover 
exponent ~ 1.78(6), which is argued to be equal to i^y, the relaxation time exponent of the DP 
universality class on a general footing. This conjecture is also confirmed in the case of the directed 
Ising (parity-conserving) class. Finally, we discuss the effect of diffusion to the IM-type models and 
suggest an argument why diffusive models with some hybrid-type reactions should belong to the DP 
class. 

PACS numbers: 64.60.Ht,05.70.Ln,89.75.Da 



I. INTRODUCTION 

The directed percolation (DP) has been studied ex- 
tensively as one of typical dynamic critical phenomenon 
far from equilibrium yy]. Nonequilibrium phase transi- 
tions of systems with a unique absorbing state are found 
to belong to the DP class if no symmetry or conserva- 
tion of the order parameter is present . Even systems 
vifith infinitely many (IM) absorbing states such as the 
pair contact process (PCP) ^ are believed to share the 
same critical behavior with the DP, but the theoretical 
understanding of the IM-type DP (DPjm) models is still 
lacking. For example, the phenomenological field theory 
introduced in Ref. [4] and elaborated in Ref. Q shows in- 
consistency with the numerical studies. In fact, the field 
theory based on the phenomenological Langevin equation 
predicts that the PCP should belong to the dynamic per- 
colation rather than the DP class [5]. Besides the uni- 
versality issue, its spreading dynamics is also not fully 
understood [1, 

One may reduce the number of absorbing states signif- 
icantly by introducing particle diffusion (PCP with diffu- 
sion or PCPD) [8]. Surprisingly, the PCPD has brought 
up serious turmoil, which could not be settled down in 
spite of extensive numerical 0, [l^, [HI and analytical 
|12l | studies. The answer seems to be one of two possi- 
bilities: The PCPD belongs to the DP class with a long 
transient, or forms a new universality class distinct from 
known universality classes to date. As an attempt to 
resolve the issue, the present authors suggested two dif- 
ferent approaches to the PCPD. 

First, introducing the dynamic perturbation which is 
implemented by the biased hopping, we showed that 



the one-dimensional PCPD with biased diffusion (driven 
PCPD or DPCPD), exhibits critical scaling distinct from 
the unbiased PCPD, but instead shares the critical be- 
havior with the two-dimensional PCPD without bias (di- 
mensional reduction) Since all known DP models 
are robust against the biased diffusion [l^, a DPCPD- 
type variant can serve as a litmus test for the PCPD scal- 
ing [l^ . Second, we studied the crossover behavior from 
the PCPD to the DP by introducing single-particle anni- 
hilation/branching reactions and showed that there are 
diverging crossover scales with the universal nontrivial 
crossover exponent 16']. These results provided another 
evidence that the PCPD is distinct from the DP. 

In this paper, we study the crossover behavior from 
the DPiM to the DP in order to understand better the 
difference between these two "equivalent" (DP) univer- 
sality classes. Actually, it would be absurd to talk about 
the crossover between two models belonging to the iden- 
tical universality class (see Sec. |lT]). However, the DPjm 
models differ from the DP models in regard to the "non- 
order" parameter at the transitions, which shows a dis- 
continuous singularity at the crossover. This singularity 
induces a well-defined and nontrivial crossover from the 
DPiM to the DP, more generally a crossover arising from a 
considerable reduction of the number of absorbing states 
(between two different DPim)- 

We find two distinct crossover behaviors depend- 
ing on the routes to reduce the number of absorb- 
ing states. Infinitesimal inclusion of an "excitatory" 
process (like single-particle branching) makes the sys- 
tem overly active, which gives rise to a discontinuous 
crossover (see Sec. lIII|) . While, the opposite "inhibitory" 
route (like single-particle annihilation) reveals a contin- 
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uous crossover with the nontrivial crossover exponent 
(/) ~ 1.78(6). We find that this crossover exponent is 
universal for various kinds of models including the PCP 
and the triplet contact process (TCP). We argue that the 
crossover exponent is not independent but equal to v\]^ , 
the relaxation time exponent of the DP universality class 
on a general footing. This conjecture is also confirmed 
in the case of the directed Ising (parity-conserving) class 
(see Sec.|lV|. 

The PCPD per se can be considered one of crossover 
models from the DPim by allowing diffusion, which re- 
duces the number of absorbing states considerably (from 
exponentially many to linearly many absorbing states 
with respect to system size). However, the crossover 
study from the PCP to the PCPD does not give any 
useful information on the DPim, because the particle dif- 
fusion makes the system more active (excitatory) and the 
crossover turns out to be discontinuous. 

Finally, we study the crossover from the diffusive reac- 
tion models to the DP. Such an example is the crossover 
from the PCPD to the DP studied in Ref. In 
Sec. |Vl the hybrid-type models with diffusion (tpl2: 
2A ^ 0, 3A -> AA) [3 are perturbed by adding a single- 
particle annihilation process {A — > 0) and its crossover to 
the DP is investigated. These hybrid-type models (where 
the branching process is of higher order than the anni- 
hilation process) are numerically known to belong to the 
DP class. As expected, the critical line emanates "lin- 
early" from the hybrid- type model point. We suggest 
an argument why these models should belong to the DP 
rather than a PCPD-type nontrivial class. 



II. CROSSOVER BETWEEN THE IDENTICAL 
UNIVERSALITY CLASS? 

This section considers a d-dimensional stochastic sys- 
tem of hard core particles with dynamics summarized 
in Table HI In Ref. [13] , it is shown that two different 
stochastic systems modeled by tilded and untilded rates 
are equivalent if transition rates satisfy the relations 

^ o" „ O" 
D + - = D + 

~ ~ T (1) 

^+i + 2 = ^(^+i + 2)' 

\ + fl + a~2iD^\ + 'q + a-2D, 

where any positive number which renders all tilded 
and untilded transition rates be nonnegative is physically 
meaningful. By equivalence is meant that all correlation 
functions of the model with tilded parameters can be de- 
duced from those of the model with untilded parameters 
and vice versa. For instance, the particle density p at 
time t of two different stochastic many body systems be- 



TABLE L d-dimensional reaction-diffusion processes of single 
species with hard core exclusion and their rates. 



Diffusion 


^ 




with rate D/d 


Pair annihilation 


AA- 


>00 


with rate \/d 


Coalescence 


AA^ 




with rate r)/ {2d) 


Coalescence 


AA^ 


%A 


with rate 77/ (2d) 


Death 


A 





with rate ^ 


Branching 




AA 


with rate g /{2d) 


Branching 


yl0 ^ 


AA 


with rate a/ (2d) 



comes [T3| 

p[D, A, i-j, cr; pq, t) = -p{D, A, fj, |, a; po, t), (2) 
M 

if the initial density of two systems has the relation po — 
fipQ. Needless to say, both po and po should lie between 
and 1. Although the equivalence is shown only for the 
one-dimensional systems in Ref. , Eq. ([1]) is generally 
true for any d-dimensional systems, which can be easily 
shown by the same technique developed in Ref. 

Now consider the branching annihilating random walks 
with one offspring (BAWl) 18] which corresponds to the 
model with ^ = in Table |T1 The parameters used in 
Ref. [1^ in one dimension are D — p/2, A — p, and 
T] = a = 1 — p with the tuning parameter p. If w = p — 1 
is very small and nonnegative, the solution of Eq. ([T]) up 
to the order of w is 

~ w ^ 

D — D = (T, a — (7 — wa, £ = wa, 

2 ' (3) 

A - A ~ -7«(ry + 2A), ij - i] ~ w{ri + 2X ~ 2(j). 

If w is sufficiently small, it is always possible to associate 
the BAWl with a stochastic process with spontaneous 
death in an equivalent way with all nonnegative rates. 
Since the BAWl in high dimensions is also known to have 
a nontrivial transition point [l^ , the following discussion 
is valid in any spatial dimension. The transition points 
for stochastic systems with small w can be always calcu- 
lated exactly from Eq. ([T|) (approximately from Eq. ([3])), 
if the transition point of the BAWl is given. 

The conclusions from the above analysis are two-fold. 
First, it is clear from Eq. ^ that the phase boundary 
(critical line) should meet the BAWl transition point lin- 
early with finite slope, as w vanishes. This implies that 
there is no additional singularity involved near w = 0, 
which is fully expected from the crossover between mod- 
els with the identical universality class. If one defines 
the crossover exponent (j) from the shape of the criti- 
cal line near w = (see Sec. IIII[) . one can say that 
0=1. Second, the critical decay of the density is given 
by Pc{t) = {1 + w)pc{t) from Eq. ([2]), which implies that 
there is no diverging crossover time scale for small w. 

Since the introduction of the spontaneous death does 
not change the structure of the absorbing phase space 
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(single absorbing state) let alone the universality class, 
the above analysis is in good harmony with the naive 
expectation as to the "crossover" between two models 
belonging to the identical class. In the next section, 
however, we will show that the substantial change of the 
absorbing phase space without affecting the universality 
class will trigger a nontrivial crossover. 



III. CROSSOVER FROM THE DPim TO THE 
DP 

Unlike the BAWl, the pair contact process (PCP) is 
the prototype of the DPjm models with exponentially 
many absorbing states. By introducing single-particle 
reactions to the PCP, the number of absorbing states 
changes drastically from infinity to one. This section 
shows that this qualitative change is reflected into the 
singular behavior of the phase boundary close to the PCP 
transition point in one dimension. 

The dynamics of the model is summarized as 

AA ^ 00, ^^^1 AAA, (4a) 

f I i^^^ 00, f I ^ AA, (4b) 
%A\ 0AJ ' ^ ' 

where < g < 1. For the PCP at w = 0, any configura- 
tion without a pair of neighboring particles (a mixture of 
isolated particles and vacant sites) is absorbing and its 
number grows exponentially with system size. The order 
parameter of the PCP is the pair density (the number 
density of AA pairs) and the particle density field is aux- 
iliary which is finite even in the absorbing phase. At 
nonzero w, an isolated particle becomes active and only 
the vacuum becomes the true absorbing state. In this 
case, the particle density is usually adopted as the order 
parameter and the pair density scales in the same way. 

Figure [1] locates the transition point of the PCP {w = 
0) at po = 0.077 0905(5) by exploiting the critical decay 
of the pair density as Pp{t) ~ with 5 to be the crit- 
ical exponent of the DP class whose accurate value can 
be found in (20l |. In numerical simulations, the system 
size is L = 2^® and the number of independent samples 
are 750, 1500, and 400 for the data in the active, criti- 
cal, and absorbing phases, respectively. The flatness of 
Pp{t)t^ over four log decades in time confirms the solid 
DP critical scaling of the PCP. 

At finite w, the model still belongs to the DP class 
irrespective of q. Unlike the PCPD to the DP crossover 
model in Ref. [Iq, however, the critical lines show two 
completely different singular behaviors, depending on the 
value of q. For large q, the activity of the system is 
enhanced by additional single-particle reaction processes 
(excitatory process) and the system becomes overly ac- 
tivated even with infinitesimal w. The critical line does 
not converge to the PCP critical point as w decreases to 




10° 10^ 10^ 10^ 10* lo'' 10*^ lo'' 
t 



FIG. 1: (Color online) Semilogarithmic plot of Pp{t)t vs t of 
the PCP near criticality with 5 = 0.1595 to be the exponent 
of the DP. Since the upper (lower) curve veers up (down), 
we estimate the critical point as po ~ 0.077 0905(5) with the 
error in the last digit by 5. 



zero (w = 0^) and shows a discontinuous jump. On the 
other hand, for small q, the system activity is suppressed 
(inhibitory process) and the system becomes more inac- 
tive. The critical line nicely converges to the PCP critical 
point and shows a continuous crossover with a nontrivial 
crossover exponent. 

First, we choose the q = 1 case as a typical excita- 
tory route of the crossover from the DPim to the DP. As 
shown in Fig. [21 the critical line approaches p ~ 0.1448 
as w approaches zero, which is by far above the criti- 
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FIG. 2: (Color online) Plots of p{t)t^ vs t ior w = 10"^ 
10~*, and w — 10~^ at q — 1 close to the critical points 
in semi-logarithmic scales. Again, 5 assumes the DP value. 
The curves in the middle correspond to p = 0.1451, 0.1448, 
and 0.1448 (from bottom to top), respectively and the value 
of p's of other two curves are ±0.0001 off from the middle 
value. As w becomes smaller, the relaxation time becomes 
larger though the critical point does not change much and 
approaches p ~ 0.1448. 



cal value of the PCP (po ~ 0.077). So there is a big 
jump of the critical line at w = 0. The discontinuity 
can be understood as follows: Consider a system with p 
slightly above the PCP critical point po and < w ^ 
where r is the relaxation time which is finite off critical- 
ity. Then the single-particle branching event {A 2A) 
with the characteristic time of occurs effectively af- 
ter the system falls into one of the PCP absorbing states 
in which the isolated particle density is finite. Since the 
branching event creates a new pair, the system is reac- 
tivated and performs the damage-spreading-type "defect 
dynamics" for some time proportional to t and again 
falls into one of the PCP absorbing states. This defect 
dynamics continues forever with the period of time . 
As the particle density is finite (and quite large) even 
in the PCP absorbing states, the time-averaged particle 
density in this iterated process should be finite in this 
region of the phase diagram. This implies that the con- 
tinuous absorbing phase transition at infinitesimal w into 
vacuum should occur way above po, which is consistent 
with our finding. Note that the discontinuity in the aux- 
iliary field (particle) density is crucial in this crossover. 

Actually, the same argument can be applied to the 
crossover from the PCP to the PCPD. We can introduce 
the diffusion rather than the single-particle reactions and 
again consider p slightly above pq. Let po denote the iso- 
lated particle density at the PCP absorbing states, then 
the characteristic length scale between isolated particles 
is 1/po- If < Dpq -C with the diffusion constant 
D, the "defect dynamics" will continue again indefinitely 
for small D. So the phase boundary in the D — p plane 
should have a discontinuity at D = 0. 

Let us turn to the crossover model with q = 0, which 
should represent a typical inhibitory route. Table [Til sum- 
marizes the critical points of the model for some w's at 
q = and the corresponding phase boundary is plotted 
in Fig. [21 Unlike the previous case, the reactive phase 
shrinks continuously with the rate of additional single- 
particle annihilation process and the phase boundary is 



TABLE II: Critical points of the model with dynamics of 
Eq. Q for some values of w's at q = 0. The numbers in 
the parentheses indicate the error of the last digits. 



w 


Pc{w) 





0.077 0905(5) 


10"^ 


0.077 002(3) 


5 X 10"^ 


0.076 885(3) 


10-* 


0.076 784(4) 


2 X 10"* 


0.076 642(2) 


3 X 10"* 


0.076 530(5) 


4 X 10-* 


0.076 432(2) 


5 X 10-* 


0.076 345(2) 


6 X 10"* 


0.076 264(1) 


10"^ 


0.075 988(4) 
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FIG. 3: (Color online) Phase boundary of the model of Eq. (|4]) 
at g = in w — p plane. Symbols locate the numerically esti- 
mated critical points. The error of the critical point is smaller 
than the symbol size. The curve shows the least-square-fit re- 
sult of the phase boundary. The absorbing (active) phase is 
above (below) the curve. Inset: the same but the vertical axis 
is Ac(io) in log-log scale. The slope corresponds to the inverse 
of the crossover exponent which is estimated as (p = 1.78(6). 

continuous. The usual analysis method can be applied 
to this case [2l[. If we define A = {po — p)/po and 
Ac{w) — (po — Pc{w))/po, the phase boundary is well 
fitted by - w^/"^ with = 0.56(2) or = 1.78(6); 
see the inset of Fig. [31 Let us assume the existence of the 
well-defined crossover scaling which is described by the 
scaling function [2l[ 

Pp{A, w; t) ^ r^TiA^i^ t, w^ii t), (5) 

where pp is the pair density and = with es- 

timated in the above. We examine whether the scaling 
function in Eq. ^ correctly describe the crossover near 
the PCP critical point. 

First, we measure the pair density for various w's at 
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FIG. 4: (Color online) Log-log plot of the scaling function 
Eq. (|6]l for the PCP crossover model using 5 — 0.1595 and 
p = 0.97. All curves are collapsed into a single curve. 
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FIG. 5: (Color online) Scaling plot of Pp{t)r vs w'^ll t with 6 = 
0.1595 and /iy = 0.97 at criticality in semi-log scales. Inset: 
Scaling function Ti. at finite w and w = 0. The asymptotic 
value of = is different from the w limit. 
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FIG. 7: (Color online) Scaling plot of pt{t)r vs w'^H with 
the same exponents in Fig. |4] for the crossover TCP model of 
Eq. ([51 in semi- log scales. As in Fig. 3] all curves are collapsed 
into a single curve. 



A = 0. From the scaling ansatz ([5]), the pair density at 
A = should collapse as 

t^Ppit)=g{w^H). (6) 

With /ijj ~ 0.97, all curves for the pair density are col- 
lapsed into a single curve as Fig. 2] shows. Next, we take 
A = Ac(w) along the critical line. Since Ac(w) ~ w^/"^ 
and = the scaling function should take the form 

t^PpiA,{w);t)^n{w^H), (7) 

where Ti{x) approaches a constant as x — > oo. In Fig. [51 
all curves at different critical points collapse well into a 
single curve. Hence we conclude that the scaling function 
Eq. ([5]) correctly describes the crossover behavior from 
the DPiM to the the DP. 




t 

FIG. 6: (Color online) Log-log plot of the particle density 
for some finite w's at A = 0. For any finite value of w, the 
particle density p{t) will go to zero, which is clearly different 
from the w = case. 



Since the models at w = and at w ^ belong to the 
same DP universality class, it is natural to ask what is 
the origin of such a nontrivial singularity near the PCP 
critical point. The inset of Fig. [H] gives a hint to this 
question, which shows that the scaling function Ti. (the 
amplitude of the critical decay) does not approach the 
PCP value as w goes to zero, i.e., it is not continuous at 
w = 0. The discontinuity in this amplitude must orig- 
inate again from the discontinuity in the auxiliary field 
(particle) density. One can see it directly from the be- 
havior of the particle density (p). Unlike the pair density, 
p can not be described by the scaling function ([5]). Con- 
sider again the case at A = and nonzero w. For any 
finite value of w in the thermodynamic limit, p{t) ap- 
proaches to zero as < cxo; see Fig. [51 On the other 
hand, the model at A = w = (the critical PCP), has 
a nonzero density of p(t) as i — > oo. In other words, the 
w — > limiting process is different from the w = model 
itself in regard to the auxiliary field density. 

To check the universality of the crossover exponent, 
we study the modified PCP with the replacement of 
2A 3A with 3^ 4A in Eq. (g]) which is the model of 
Eq. ([9]) with no diffusion {D — 0). This model also has 
infinitely many absorbing states and belongs to the DPim 
class. By introducing single-particle reactions {w ^ 0), 
the same crossover behavior is found as the above (data 
not shown). 

We also study more general crossover behavior from 
one DPim to another DPim with the considerably re- 
duced number of absorbing states. To be specific, we con- 
sider the triplet contact process (TCP) and its crossover 
model by introducing the 2A — ^ A process without spon- 
taneous death. The TCP with pair dynamics is defined 
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AAA ^ 000, ^^^^1 (^-rt/^ AAAA, 



' 0AAAJ 



(8) 



AA0 



i«/2 



0AA 00A. 



The above model has infinitely many absorbing states, 
but with nonzero w the number of absorbing states is 
greatly reduced. At w = 0, there is again a jump in 
the auxiliary field density (here, the pair density). We 
found that the critical point for the TCP at w = is 
Pc = 0.036 865(5), exploiting the DP critical scaling (data 
not shown). Figure [7] shows the scaling plot of the triplet 
density pt in the same way as in Fig. U) We also mea- 
sured the crossover exponent from the phase boundary 
and found the same exponent (data not shown). 

Hence we conclude that there is the well-defined and 
universal crossover scaling from the DPjm to the DP 
which is mediated by the significant reduction of the 
number of absorbing states. The discontinuity in the aux- 
iliary field plays a crucial role in this nontrivial crossover. 



IV. CONJECTURE ON THE CROSSOVER 
EXPONENT (j) 

The crossover exponent from the DPjm to the DP is 
estimated as (/) = 1.78(6). Since this crossover occurs 
between the same universality class, we are suspicious 
that (f) may not be independent but related to the well- 
known DP critical exponents. Actually, we argue that 
the crossover exponent is given by the DP relaxation 
time exponent: (p — v\\ ~ 1.733 (or /iy — 1) which is 
compatible with the numerical estimation within error. 
The reason is as follows: Take the model of Eq. @ at 
(7 = 0. The critical line should be determined by the 
competition between the single-particle annihilation pro- 
cess (A — > 0) parameterized by w and the multi-particle 
(pair) reaction process (2A or 3A) parameterized 
by A (po —p)- We expect both events should appear 
at the same time scale along the critical line to balance 
off each other. Since the single-particle (auxiliary field) 
density is finite at the PCP (DPjm) point, the time scale 
for A ^ should be simply proportional to w^^ . The 
time scale for the pair reaction process should be given by 
the relaxation time scale r ~ A""^!! . Consequently, the 
critical line is determined as Ac ^ {pa ~ Pc{w)) ~ u'^^'^ii , 
which yields (j) — and equivalently — 1. 

Considering the crossover between the DP models with 
a unique absorbing state, the time scale for the process 
parameterized by w is proportional to like the other 
competing process because the auxiliary field (particle) 
density is also vanishing critically as w decreases to zero. 
In this case, we get 0=1 and equivalently = 
which is consistent with our result in Sec. HTl 

Since our argument for the crossover exponent is gener- 
ally applicable to any universality class, we can check its 



validity through studying the similar type crossover be- 
tween the directed Ising (DI) class models |22|] . Consider 
a one-dimensional system with two species, say A and B. 
Between the same species, hard core exclusion is applied, 
but different species can reside at the same site. The dy- 
namic rules are as follows: The dynamics always starts 
with an A particle. A randomly chosen A particle can 
hop to one of nearest neighbors with probability p. If two 
A particles meet at the same site by hopping, both parti- 
cles are removed with probability i. If this annihilation 
attempt fails, the particle goes back to the original site. 
With probability 1 ~ p, a B particle is generated at the 
same site occupied by the chosen A particle. If that site is 
already occupied by another B particle, the two B parti- 
cles transmute to two A particles which will be placed at 
two nearest neighbor sites. If any of transmuted A par- 
ticles is placed at the site already occupied by another A 
particle, both particles are annihilated immediately. In 
summary, 2A ^ 0, A ^ A -I- i?, and 2B 2A processes 
are allowed with A particle diffusion. 

Since B particles are not allowed to hop, the system is 
inactive without an A particle but only with B particles. 
The number of the absorbing states grows exponentially 
with system size and the auxiliary field (B particle) den- 
sity is finite at the absorbing transition. Besides, the 
number of A particles is conserved modulo 2, which is 
the characteristic of the DI (or parity-conserving) class. 

As in Sec. IIIIl we study the crossover by introduc- 
ing spontaneous annihilation of B particles (inhibitory 
route) with rate w. As expected, we find the DI criti- 
cal scaling for both w = and w cases. We locate 
the critical line by exploiting the known DI critical expo- 
nents [22I, which is summarized in Table Hill for some w's. 
Since the critical exponent i^n of the DI class (~ 3.25) is 
much larger than that of the DP (~ 1.733), the accu- 
racy of the critical points in Table IIIIl is worse than that 
for the DP cases. From these data, one can estimate 
the crossover exponent 1/(/)di as 0.34(4) which should be 
compared with l/i^n of the DI class (~ 0.31). Hence we 
conclude that our argument for (j) ~ v\\ also applies to 
the crossover from the IM-type DI to the DI models. 



TABLE III: Critical points of the DI crossover model for some 
values of ui's (see the text). The numbers in the parentheses 
indicate the error of the last digits. 



w 


10-* 
3 X 10"* 

10"^ 
3 X 10"^ 



Pe(w) 



0.4518(1) 
0.4443(2) 
0.4405(2) 
0.4353(1) 
0.4276(1) 
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FIG. 8: (Color online) Semi-log plot of p{t)t^ vs t with 
5 = 0.1595 at pc = 0.058 427, 0.058 723, 0.058 966, and 
0.059 031 for w = 10"^ 5 x 10"*, 10""*, and (from left 
to right), respectively. There is no diverging time scale as w 
approaches zero. Inset: Phase boundary near w = 0. Each 
symbol corresponds to the critical point used in the main 
figure. The straight line is drawn between two consecutive 
points, independently. 

V. DIFFUSION EFFECT 

This section studies how the crossover scaling is af- 
fected if particles are allowed to hop in models consid- 
ered in Sec. IIIII With single-particle diffusion, the PCP 
becomes the PCPD where the particle (auxiliary field) 
density as well as the pair density vanishes at criticality 
even without any single-particle reaction process. The 
crossover from the PCPD to the DP caused by including 
single-particle reactions has been studied previously by 
the present authors [IB] where the value of the crossover 
exponent is reported as l/<j) = 0.58(3). This value is 
quite close to that obtained for the crossover from the 
DPiM to the DP in Sec. EH This similarity may mislead 
one to jump to the wrong conclusion that the critical na- 
ture of the PCP and the PCPD is equivalent. However, 
one should remember that the origin of the nontrivial 
crossover from the PCP to the DP lies in the finiteness 
of the auxiliary field density at the PCP critical point, 
while it vanishes at the PCPD critical point. So, if the 
PCPD belongs to the DP class, then one should expect a 
trivial crossover with (j) — \. Our finding of the nontrivial 
crossover in [l^ implies that the PCPD class is distinct 
from the DP class. Hence it is likely that the similarity 
of two crossover exponents be a mere coincidence. 

The distinction of the PCPD from the DP can also be 
evidenced by the study of the crossover model with the 
hybrid-type reaction dynamics 2A — > and iA AA. 
Without diffusion, this model belongs to the DPim and 
its crossover to the DP was studied in Sec. IIIII Even 
if particles are allowed to diffuse, this model (called as 
tpl2) is numerically known to belong to the DP class @. 

It will be illuminating to see how the diffusion in the 
tpl2 can change the crossover behavior to the DP. The 



dynamics of the model is summarized as 

(9) 

M 00, A^^U, 

9Aj 

where D = (1 ~ w)/2. The model at w = is the 
tpl2. Our numerical results in Fig. [8] show the typical 
"crossover" behavior between the identical universality 
class discussed in Sec. [TTl The critical line converges to 
the u; = point linearly {(j) — 1) and there is no diverg- 
ing time scale as w becomes smaller with almost perfect 
collapse of all critical density decay curves near w k, {). 
Our crossover study provides another strong numerical 
evidence that the tpl2 belongs to the DP class. 

If the PCPD does not belong to the DP class, why 
should the tpl2 belong to the DP class? The PCPD and 
the tpl2 seem quite similar in the sense of multi-particle 
nature in reaction dynamics and also the absorbing space 
structure (vacuum and a single-particle state). However, 
they are quite different in the role of diffusing isolated 
particles. See Fig. [S] for the space-time configurations for 
the tpl2 and the PCPD at criticality, starting from the 
low density initial condition without pairs. The diffusing 
isolated particles of the tpl2 cannot increase the num- 
ber of particles in most cases, because 2A — *■ dynamics 
dominates over 3A — > 4A dynamics: Pairs generated by 
collisions of two isolated particles evaporate before greet- 
ing another isolated particle to become "active" triplets. 
Consequently there is effectively no feedback mechanism 
from isolated particles to increase the particle density or 
make the system more active. Therefore the region of iso- 
lated particles can be regarded as absorbing like in the 
PCP model. This case may correspond to the no-feedback 
point (r = 0) for the generahzed PCPD (GPCPD) stud- 
ied in [IQ] , which is the DP fixed point. This argument 
can be generalized to systems with hybrid-type reaction 
dynamics mA (m + k)A and nA —>■ {n — l)A with 



X X 




FIG. 9: (Color online) Space-time configurations of (a) the 
tpl2 and (b) the PCPD at criticality. The initial density is 
(a) ^ for the tpl2 and (b) for the PCPD with size L = 2^°. 
Isolated particles (particles in a cluster) are designated by the 
blue (red) color. 
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m > n and k,l > 0, which are numerically shown to 
belong to the DP class 0. 

The isolated particles of the PCPD, however, cannot 
be regarded as absorbing as Fig. [5] shows; the isolated 
particles may affect the critical spreading actively, be- 
cause both dynamics of 2A and 2A — + 3A compete 
each other and consequently there is an effective feedback 
mechanism from diffusing particles to make the system 
active. This corresponds to the GPCPD with long-term 
memory effects at r ^ [13]. 

VI. SUMMARY AND CONCLUSION 

In summary, we studied the crossover from the model 
belonging to the directed percolation (DP) class with in- 
finitely many absorbing states (DPjm) to the DP model 
by reducing the number of absorbing states significantly. 
The crossover is found to be well described by the usual 
crossover scaling function for the order parameter. The 
crossover exponent (j) is argued to be related to one of 
the DP critical exponent, i.e., = which is further 
evidenced by the similar crossover model belonging to 



the directed Ising class. The origin of the diverging scale 
and the nontrivial crossover comes from the discontinuity 
of the auxiliary field density at the DPjm critical point. 
Our study for the first time presents the existence of the 
nontrivial scaling in the DPjm, which is compared with 
the study on the spreading exponents. 

We also studied how the crossover behavior from the 
DPiM to the DP is affected by particle diffusion. The 
crossover from the pair contact process with diffusion 
(PCPD) to the DP studied in Ref. [H] is well classi- 
fied by the nontrivial crossover exponent. On the other 
hand, the tpl2 which is known to belong to the DP class 
is characterized by the trivial "crossover" between the 
identical class detailed in Sec. |TT1 This provides an ad- 
ditional evidence supporting that the PCPD is distinct 
from the DP. In addition, we suggest an argument based 
on the role of diffusing isolated particles why the tpl2 
should belong to the DP class, but the PCPD does not 
need to be. 

It will be a challenging problem to see if the crossover 
scaling from the DPim to the DP can be anticipated in 
the framework of the field theory 
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